Landau— Zener interferometry for qubits 
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One may probe coherence of a qubit by periodically sweeping its control parameter. The qubit 
is then excited by the Landau-Zener (LZ) mechanism. The interference between multiple LZ tran- 
sitions leads to an oscillatory dependence of the energy absorption rate on the sweeping amplitude 
and on the period. This interference pattern allows to determine the decoherence time of the qubit. 
We introduce a simple phenomenological model describing this "interferometer" , and find the form 
of the interference pattern. 



During the last few years, a number of proposals 
for constructing quantum bits (qubits) from mesoscopic 
Josephson junctions have appeared Q H, Q an d 
first experimental results in this direction have been re- 
ported 0, IS S B 0> El EH]- Large part of these 
qubits are actually different physical realizations of an 
externally controllable quantum double- well system, with 
nearly equal depths E\ t 2 of both wells \E\ — E^\ <§C luq <C 
Ei (here loq is the oscillation frequency within a sin- 
gle well), and with the inter- well tunneling amplitude 
A ~ \Ei — E%\. The above conditions ensure that higher 
eigenstates of the system are separated from the nearly 
degenerate doublet by a large gap (compared to A), and 
the probability of their excitation can be neglected. The 
energy difference \E\ — E?\ is controlled by an external 
time-dependent parameter x(t), which is either voltage 
for the SET-bascd "charge" qubit 0, 0, or magnetic 
flux through the Josephson junction loop for the "phase" 
qubit y>y,|a,l3- A review of recent results for both types 
of superconductive qubits can be found in [12j . Quantum 
manipulations with qubits involve varying in time both 
x(t) and A(t) (as well as more complicated two-qubit 
manipulations). Since the overall scale of possible E\^, 
variations as function of control parameter X is very large 
compared to relevant values of A (e.g. it was more than 
100 times larger in the design of Ref. |2(), fluctuations of 
x(t) are expected to be one of the most important sources 
of dephasing in such qubits Q . Indeed, the experimental 
data of Ref. seem to confirm these expectations. 

Thus the first problem to be addressed in the develop- 
ment of this type of qubits is to find a convenient probe 
which tests whether the device undergoes coherent evo- 
lution. Resonant absorption method was used in experi- 
ments of Refs. 0,0], whereas Nakamura et al have used 
time-domain manipulations Q. It may be preferable to 
employ simpler methods to measure decoherence time of 
a qubit, without super- high-frequency (in the GHz range) 
manipulations. One such method (based on the measure- 
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FIG. 1: Eigenvalues e of the Hamiltonian Q depend 
on the control parameter x. When the control parameter 
passes the Landau-Zener point (x—0), the qubit may be non- 
adiabatically excited with the probability v. The upper state 
may also decay into the lower one, dissipating the energy into 
the environment. The decay probability is w. 



ment of static nonlinear Andreev conductance) was pro- 
posed in Ref. [13| for the specific case of superconducting 
phase qubit like those proposed in Refs. Another 
low- frequency probe is the observation of Ramsey fringes, 
which was performed for superconducting qubits in 0,0] • 
In the present paper we propose and analyze a different 
method to determine the decoherence time of a qubit by 
a low-frequency non- resonant measurement. 

The idea is to employ the Landau-Zener (LZ) non- 
adiabatic tunneling processes. If the control parameter 
x(t) changes in such a way that the "collision region" 
with A ~ \E\ — £2 1 is traversed, the system may become 
excited from the lower to the upper level. If the pa- 
rameter x{t) is changed periodically, successive Landau- 
Zener tunneling events interfere, and such an interference 
allows to estimate the decoherence in the qubit. The 
interference picture is easiest to describe in the setup 
with the amplitude of the parameter sweep x large com- 
pared to the width of the "transition region" (Fig. ^| . If 
the dephasing time ideph is larger than the period of the 
parameter sweep to, the interference between successive 
Landau-Zener transitions has an oscillatory dependence 
upon xq and to- The interference may be observed from 
the rate of energy dissipation which is proportional to 
the occupation of the upper level. The precise form and 
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the strength of the interference pattern is determined by 
the interplay of the coherent Landau-Zener transitions 
and of the spontaneous incoherent decay E 2 — ► E\. Per- 
forming interference experiments at a constant offset x Q g 
in the oscillations of the parameter x(t) and at differ- 
ent values of xq and t allows to determine the rates of 
dephasing and of inelastic relaxation in the qubit. 

Specifically, we model the qubit by the two-level 
Hamiltonian 



H 



x(t) + i(t) 



- x (t) - at) 



(i) 



In this paper, we consider only one mode of qubit opera- 
tion, namely varying the diagonal controlling parameter 
x{t) while keeping the gap A fixed. We also assume that 
the main source of decoherence are thermal and quantum 
fluctuations of the control parameter around its intended 
value. Such fluctuations are described by the operator 
£(i). Phase fluctuations of A may also be incorporated 
into £(i) by an appropriate gauge transformation. 

An example of the physical system leading to the 
Hamiltonian Q is the superconducting phase qubit 
0, 0, 0| . The control parameter x(t) in this design is 
the magnetic flux through the qubit loop. The diago- 
nal coupling to the environment is realized via coupling 
to the magnetic flux (including the fluctuations of the 
external electromagnetic field 0, 0] and coupling to nu- 
clear magnetic spins ^^)- We also remark here that 
the diagonal coupling to the environment in the Hamil- 
tonian JQ) does not represent the most general form of 
coupling. In general, off-diagonal coupling may also be 
present, which would lead to the inelastic decay of the 
qubit states away from the Landau-Zener transition re- 
gion. However, we assume that the off-diagonal coupling 
is already sufficiently suppressed: the qubit can preserve 
its occupation-number information in the "idle" regime 
(away from the Landau-Zener transition region), and the 
decoherence is determined by the diagonal coupling chan- 
nel. 

The quantum variable £(i) describes coupling to the 
collective degree of freedom of the external reservoir. It 
is a quantum variable corresponding to a collective de- 
gree of freedom of the reservoir. The usual model for 
the reservoir is an ensemble of harmonic oscillators [jfij. 
Our treatment will be phenomenological and not involv- 
ing the microscopic properties of the reservoir, therefore 
we do not explicitly include the reservoir in the Hamilto- 
nian Q. However, we check the validity of our approach 
by comparing it to the microscopic calculation for the 
oscillator bath model (see Appendix). 

Depending on the experimental conditions, the tem- 
perature of the reservoir may be either lower or higher 
than the gap A. We first consider the case of the reser- 
voir temperature much smaller than the gap A, and later 
explain how the results are modified at higher tempera- 
ture. Independently of the relation to the gap A, we 



assume that the reservoir temperature is always higher 
than the sweep frequency t$ : this is necessary for our 
treatment of dephasing as a Gaussian noise and for our 
assumption of independent dephasing processes on dif- 
ferent half-periods of the parameter sweep. 

The effect of the coupling to £(£) is twofold. In the 
transition region (x(t) ~ A), this coupling has non- 
vanishing matrix elements between the two adiabatic lev- 
els, and therefore leads to inelastic transitions between 
the levels. In the limit of the reservoir temperature 
much lower than A, the transitions occur mostly from 
the upper to the lower level, thus attenuating the tran- 
sition probability |l4| . Away from the transition region 
(|x(i)| ^> A), the Hamiltonian is almost diagonal, 
and the effect of £(i) is dephasing. 

The control parameter x(t) is swept periodically, with 
the amplitude xq and with the period t$: 



. 2irt 

x(t) = xo sin h x g 

to 



(2) 



For simplicity, we first consider in detail the case of zero 
offset x s = 0, and then discuss the interference pattern 
at arbitrary x g . 

Each time the control parameter passes the Landau- 
Zener point x(t) = 0, the Landau-Zener tunneling oc- 
curs. This tunneling is a quantum-mechanical process 
sensitive to the relative phase of the two states. There- 
fore, the energy absorption per period depends on the 
phase tp n picked up far from LZ point. The latter phase 
is determined by both the sweep amplitude and the fre- 
quency (assuming x Q s = 0): 



to/2 

f n = 2 J (x(t) + £(£)) dt = ip+5ipn , <p 





2x t 



(3) 



Here ip is the average phase picked up per half-period, 
and dtp n are its fluctuations. If we assume that the cor- 
relation time of £(t) is much shorter than to (which is 
equivalent to assuming that the reservoir temperature 
is much higher than t$ ), the probability distributions 
of Stp n are Gaussian and uncorrelated on different half- 
periods (labeled by the integer n). 

The amplitude xo is assumed to be much larger then 
the level-crossing region (xq ^> A), and the period to 
should be sufficiently large, so that the Landau- Zener 
transition probability |2l| is small: 



exp 



A 2 t 
2x 



< 1 



(4) 



Also, xo should not be too large so that the Hamiltonian 
JU would still adequately describe the system. (For the 
superconducting phase qubit it implies that the ampli- 
tude of the flux modulation should be small compared to 
superconducting flux quantum.) 
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We describe the system evolution in terms of the two- 
level density matrix p. The evolution per one half-period 
of the parameter variation is given by the master 
equation which includes the three effects: the coher- 
ent Landau-Zener transitions, the inelastic decay of the 
qubit, and the phase picked up during the system evolu- 
tion away from the level-crossing region. 

We separate the decoherence effect into the two parts: 
the dephasing away from the transition region and the 
inelastic relaxation in the transition region. This separa- 
tion is possible if the transition region is narrow: A <C xq. 
Instead of microscopically deriving the relevant couplings 
(see e. g. Ref. we include them phenomenolog- 

ically as independent parameters in the master equation 
on the two-level density matrix. Both types of deco- 
herence are assumed to be small. More precisely, we 
describe the decoherence by the two dimensionless pa- 
rameters: the average phase fluctuation per one sweep, 
u = (Scpn), and the probability of inelastic decay per 
one crossing of the transition region w. Our parameters 
w and u are proportional to the longitudinal and trans- 
verse relaxation rates, respectively (defined as r rc i ax and 
1 d, m Refs. EG3). These relaxation rates have contri- 
butions from different decoherence channels, their calcu- 
lation from microscopic principles is difficult, and a priori 
there is no universal relation between them. Therefore 
instead of computing them microscopically, we introduce 
them phcnomcnologically as two independent parameters 
u, w <C 1. 

The three small parameters u, w, and v depend, in 
principle, on the period and amplitude of the oscillations 
of x(t). We shall return to this dependence in the end of 
the paper. 

The Landau-Zener transition, in the absence of inelas- 
tic events, is described by the unitary rotation of the 
density matrix: 

p i-> SpSfl , S = ( r t J , (5) 

where t and r are the amplitudes of the transition and of 
staying at the same adiabatic level, respectively. Ampli- 
tudes r and t satisfy unitarity condition |r| 2 + \t\ 2 = 1. 
The magnitude of t determines the LZ transition proba- 
bility gj: \t\ 2 = v. 

Let us assume first that the reservoir temperature is 
much lower than the gap A. Then the inelastic relax- 
ation processes may be described as attenuation of the 
amplitude of being in the upper energy state (the case of 
higher reservoir temperatures will be discussed later, see 
Eqs. (EH and iJUJl below): 

p i-> V w pV w + wp m p 2 2 , (6) 

where 




The equilibrium density matrix p(°) describes the system 
in the lower energy state. The quantity P22 is the diag- 
onal element of the density matrix corresponding to the 
upper state. The first term in Eq. © describes the decay 
of the upper state. It attenuates the amplitudes of being 
in the upper energy state by y/1 — w. Thus, the proba- 
bility that the upper state will not decay is 1 — w, and 
this allows to identify w as the inelastic decay probabil- 
ity. The second term in Eq. © describes the probability 
flow into the lower state due to the inelastic decay. This 
term is diagonal since the decay is assumed to be inco- 
herent. Note that the Eq. JBJ preserves the trace of the 
density matrix. 

A complete solution of the dissipative dynamics re- 
quires simultaneously taking into account the Landau- 
Zener processes (JSJ) and the dissipative processes © sim- 
ilarly to the treatment in Refs. [l3.ll5| . which we do not 
attempt here. Instead, note that if the transition rates 
are small (v, W <C 1), the elastic and inelastic processes 
may be considered independently. Thus, we simply com- 
bine Eqs. 10 and ©: 

p 1 ► SV wP V w S^ + w P 22P {a) ■ (8) 

The matrices V w and S do not commute, but the leading 
terms in w and v do not depend on the order of multi- 
plication. We ordered S and V w in Eq. JSJ) so that the 
trace of the density matrix is conserved. 

Finally, the phase picked up far from the Landau- 
Zener point produces the relative phase rotation of the 
upper and lower states: 

p i-> $„p$+ , $„ = cxp (iip n a z /2) , (9) 

where a z is the Pauli matrix. 

In this way, the evolution of the density matrix per one 
sweep is described by the master equation 

p h- (sV w pV w & + W P 22P {0) ) Sj, (10) 
We parameterize the density matrix as 

P=\ (ao +aa) . (11) 

The scalar part a Q — I remains constant, as required 
by the normalization of the density matrix. Then the 
dynamics is described by an equation for the polarization 
vector a. The phase factor <E>„ in rotates the vector 
a by the angle <p about z-axis. Similarly, the scattering 
matrix S in (JjjJ rotates the vector a around some axis in 
xy-plane. One may redefine the phases of the upper and 
lower states so that S describes the rotation about x axis. 
(This transformation shifts all phases ip n by a constant.) 
Thus the polarization vector a evolves per one sweep as 

a„+i = Q„a„ + wi , (12) 
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where Q n is 



Q n = R z (8 Vn )R z {ip)R x {6)Av 



(13) 



Here R z and R x are the rotation operators about z- and 
x-axes, and A w describes the attenuation: 

cos ip — sin ip 
R z ((f) — I s'mcp cos ip 
1 



R x (6) 



1 

cos 9 — sin ( 
sin 9 cos 9 



(14) 







Vl^w 

_i„ = ! o \ 1 w o 

1- w 

The parameter 9 is related to the Landau-Zener transi- 
tion probability by v = sin 2 9/2. Note that the expres- 
sion ((13(1 is correct only to the leading orders in the small 
parameters w and v and should be treated as such. 

Eq. I(12|) must be solved for a stationary solution with 
fluctuating Sip n . Since the phase fluctuations 5(p n are 
assumed to be uncorrelated, averaging over these fluctu- 
ations amounts to averaging the evolution operator Q n . 
After averaging R z (5ip n ) 



e- u / 2 



Rz(S^n) = \ e-"/ 2 
1 



(15) 



we arrive at the equation on the stationary solution a: 

a-Qa+K-z, Q = Rjfrp~n)R z (tp)R x (0)A w . (16) 

Solving 1)16(1 for a, we find the population of the upper 
level P+: 



l-a. 



vlu + w) 



\w{u + w) 2 + 2v(u + w) + 4w sin 2 ^ 

where we have kept only the leading terms in the small 
parameters w, u, and v. 

This equation describes Lorentzian peaks positioned at 
ip = 2ml. The peaks are sharp if 

w > uv, (18) 

in which case the width of the peaks Sip is given by 



2 (u + w) 2 2v{u + w) 
W) = : + 



(19) 



In the case of an arbitrary non-zero offset x g super- 
imposed onto the periodic variation of the parameter 
the two half-periods of the parameter sweep are no longer 
equivalent. The phase differences gained on odd and 
even half-periods differ by the corresponding phase offset 



ipoS = 2t Q x oS /-K: ip n — ip + 5ip n ± ip oS , with the plus 
and minus sign for even/odd half-periods respectively. 
As a consequence, the period of the master equation 1 )10(1 
doubles, as it now includes two half-periods of the pa- 
rameter sweep. In the interference pattern this produces 
secondary interference peaks at tp = n + 2nn. The rela- 
tive intensities of the two peaks depend on the offset ip Q e, 
with the two intensities equal at tp g — tt/2 + irn, and 
with one of the two peaks disappearing at ip g — 7m. 

A tedious but straightforward calculation results in the 
following extension of the formula ((17|) to the case of 
arbitrary ip Q g: 



P+{<P) = 
where 



v(u + w) [l + cosp C0Sip o g + g(tt + w) 2 



D(ip, ipos) — — w(u + w) 2 + 2w sin 2 <p 



(20) 
(21) 



+2v(u + w) 



1 + cos <p cos ip ff + — (u + w) z 



This expression is again valid only to the leading or- 
ders in the small parameters u, v, and w (and coincides 
with Eq. 1 (17(1 for ip s = only in this limit). The terms 
i(it + w) 2 in the numerator and in the denominator are 
relevant only near the points 1 + cos^> cos(p Q s = 0; away 
from these points, the terms \{u+ w) 2 are beyond the 
precision of Eq. l(2U|) and should be disregarded. Several 
examples of interference curves at different values of ip s 
are plotted in Fig. 2. Provided the condition of strong in- 
terference ((18(1 is satisfied, the height of secondary peaks 
in l(2T)|) become equal to the background at 



<PaB 



u + w 



(22) 



Note that the secondary peak is much narrower than the 
primary one as long as its height is small: the width of 
small peaks is determined solely by the strength of deco- 
herence processes u + w (only the first term in Eq. 1(19(1 ). 
whereas the width of the high primary peak involves the 
Landau-Zener amplitude v. 

So far our discussion assumed the reservoir tempera- 
ture T rcs much lower than the gap A. Taking into account 
a finite reservoir temperature, the inelastic processes in 
(|B*|) must include not only transition from the upper level 
to the lower one, but also the reverse transitions from 
the lower level to the upper one (absorbing energy from 
the reservoir). The single transition probability w should 
then be replaced by the two probabilities w± and u>2. 
Eq. JBJ is replaced by 



P i-> V w pV w + wip^ P22 + w 2 p 2 Q) pn 



(23) 



where 



V w = 









V 1 - w i 
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Pi 0) = 



1 




(0) 




1 



(24) 



We give the expressions for the transition probabilities w± 
and u>2 in terms of the environment spectral function in 
the appendix. Repeating the same derivation as before, 
we arrive to the equation 



l SL n + (wi - w 2 )i 



(25) 



replacing Eq. I|12|) , with Q n given by the same expressions 
l(T3j l and except that now 



A,, 



^(i- Wl ){i- W2 ) 

y/(i- Wl )(i-w 2 ) 

1— W\— W2 



(26) 

For small u>i and u>2, this expression for A w is equiva- 
lent to introducing the effective decay probability w = 
wi + W2- Then the solutions may be obtained from our 
previous low-temperature results by a simple rescaling 
a i—* a(uii —W2) I {w\ +W2). In terms of the average occu- 
pation number of the upper level P+(ip), this translates 
to 



P+(<p) = P+(<P,w = Wi +w 2 ) 



W 2 



Wi - W 2 



(27) 



where f +(<£>, w) in the right-hand side is given by 
Eqs. I|17|) or (|20|) . In other words, at finite reservoir tem- 
peratures, the interference pattern is simply rescaled by 
the factor (wi — W2)/(wi + W2). At low reservoir tem- 
peratures, T res <C A, the ratio w 2 /wi becomes exponen- 
tially small: W2/W1 ~ exp(— 2A/T ros ). At high reservoir 
temperatures T ros 3> A, the probabilities wi and w 2 are 
close to each other, (w\ —W2)/(wi +W2) ~ A/T rcs , which 
accordingly decreases the amplitude (but not the sharp- 
ness) of the interference pattern P+(<^j). 

Experimentally, it may be possible to measure the en- 
ergy absorption which is proportional to the population 
of the upper level P+. By observing the appearance of 
the secondary peaks at varying tp g, it should be pos- 
sible from (|22|) to determine the combined decoherence 
rate u + w. This is precisely the quantity which defines 
the quality of the qubit. The condition QlSj l should be 
fulfilled in order resolve well interference picture. Es- 
timating w ~ r rc i ax £ A/2;o and u ~ T^to: (examples of 
estimates for longitudinal and transverse relaxation rates 
Ereiajjand r<£ for superconductive qubits can be found in 
Uall^ ), and using 0, one finds the condition 



Trelax ^ I" \Xq . A to , 

> -To- exp(- 



A 



A 2 



2x 



(28) 



The condition (|28|l should be fulfilled together with in- 
equalities u,w <C 1. All these conditions together are 
compatible for low enough dephasing rate T^; taking for 
the sake of estimate T^/A = 1CP 3 and r ro i ax < (cf. 
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FIG. 2: Interference pattern P+(<p) as given by Eqs. 
and <2 U for different values of the offset ip s- The non- 
adiabatic excitation probability v = 10~ 3 , the decay prob- 
ability w = 10~ 3 , the dephasing factor u — 0.1. The four 
curves (a)-(d) correspond to ip a g = 0, vr/32, tt/8, and n/2 
respectively. 



Ref. .16] ), we find broad interval of allowed xo and 
to. Experimentally, values of T^/A ~ 10 -2 were mea- 
sured in the first superconductive qubits 0, |(J , whereas 
much smaller normalized dephasing rate of order 1CP 4 
was achieved for a non-quasiclassical device studied in 
Ref. 0. 

It may be useful to perform measurements at different 
values of the amplitude xo and period to of the param- 
eter sweep. Both Landau-Zener transition probability v 
and the inelastic decay probability w depend only on the 
velocity at the transition point. If xo and to are changed 
simultaneously so that Xo/to is kept constant, v and w 
should also remain constant. At the same time, for short- 
range correlations of £(t), the dephasing u scales linearly 
with to', u — toT^. Under these assumptions, from mea- 
surements at different xo and to it may be possible to 
determine the dephasing rate T^, and the transition prob- 
abilities v, and w. 

We are grateful to J.E.Mooij for a very useful dis- 
cussion. This research was supported by NSF grant 
PHY99-07949, NWO-Russia collaboration grant, Swiss 
NF, RFBR grant # 01-02-17759, by Russian ministry of 
science, and by the program "Quantum Macrophysics" 
of Russian Academy of Science. 



APPENDIX: THE MICROSCOPIC DERIVATION 
OF THE EQUATION FOR THE DENSITY 
MATRIX 

To establish connection between the microscopic 
Hamiltonian (JIJ and the phenomenological Eqs. (JSJl and 
(|23|l . we compute the density matrix directly. We con- 
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sider the evolution of the qubit during one sweep, treat- 
ing the coupling to the environment perturbatively, find 
correction to the density matrix p, and compare it with 
the expansion of Eq. (|23l) in small W\ and w^. First, 
we rewrite the Hamiltonian Q in the basis of adiabatic 
states: 

H = H (t) + V(t) , (29) 
where the unperturbed Hamiltonian Ho(t) is diagonal, 
H = o z t t ; e t = ^/x 2 (t) + A 2 , (30) 



x(t) 



and the perturbation V(t) is given by 

V(t) = £{t)(cosd t a z + sm6 t a x ) ; 9 t = tan" 1 

X\i 

'(31) 

We neglect the probability of the LZ transition, and con- 
sider only the transitions due to the coupling to the en- 
vironment V(t). In doing so, we assume that the char- 
acteristic energies involved are or order of et ~ A, and 
the main i-dependcnce of the perturbation V(t) is due 
to fast fluctuation of £(t). Since 9 t changes essentially 
only on the large time scale Ato/xQ, we will treat it as a 
slow variable. 

To compute the evolution of the density matrix p under 
the Hamiltonian Q, one may use Liouville equation in 
Heisenberg representation 



p{t) = i £W(t),p(t) , 



(32) 



where 



U(t) = sin0 t e^^W-2^x(') ; x (t) = 2 / e t dt , 



(33) 

where x(X) is the phase difference between the two adi- 
abatic states. The perturbation theory with respect to 
V(t) gives, to the second order, 



p(t) = p (t)+t / dt! £(ti)#(*i),A>(*) 



t u 
dt\ / dt 2 



atiWih), n(t 2 )u(t 2 ),p Q (34) 



-oo — oo 



where po is the (time- independent) density matrix in zero 
approximation. Now, we average Eq. I|34|) over fluctua- 
tions of The first order term vanishes, and one has 
for 8 pit) = pit) - pa 



Spit) 



dti J dt 2 [Ufa), U(t 2 ) Q{h - t 2 W5) 

— OO 

PoUfe) Qih - ti)] • 



Here Q(t) = (£(*)£(0)) is the correlation function of the 
environment. Note that since is quantum variable, 
Q(t) 7^ Qi~t). Rewriting the commutator in Eq. 
one finds 

t oo 

Sp n {t) = -2Re J dt x J dr Q(t) sin9 t sin6» t _ T x (36) 

-oo 

X (p ue *x(ti)-*x(ti-T) _ p22e ix(ti-r)-ix(ti)^ ( 

t OO 

Sp 12 (t) = - J dh JdT[Q{r) +Q(-t)} sin 9 t sin 6 t J£7) 

-OC 



-oo 



x I j0l2e l x(ti-r)- l x(ti) _ p 2ie ix(ti)+ix(h-T) 



and also Sp 22 (t) — —Spu(t), 8p 2 \{t) = Sp* 2 (t). The dom- 
inant contribution to the integral over r in Eq. (j36[) comes 
from the region r ~ e^ 1 , and one can use an approxima- 
tion x(ii) — x(ti — r) w 2t£ 4i . For the oscillator bath 

H cnv = ">i4ai » !(*) = E TiKe^* + a+e"^*) 

i i 

(38) 

integrals in Eq. H36J) can be expressed in terms of the 
environment spectral function \2(\ 



(39) 



as 



oo oo 

y Q(T)e lbJT dT = 2i J nj(n)dn 



n(Q) 



-fi-iO 



(40) 



where n(fi) = (exp(f2/T res ) — l) -1 is Bose-Einstein dis- 
tribution function, and T rcs is the reservoir temperature. 

After a straightforward calculation, one arrives to the 
correction to the density matrix at the end of one sweep: 

5pn(oo) = -pnw 2 +p 22 wi (41) 

<5pi 2 (oo) = -- (wi + w 2 )p 12 + i$pn . (42) 



Here 



oc 

wx = At: J [1 + n(2e t )] J(e t ) sin 2 t e t dt (43) 

— OO 
OO 

w 2 = 4tt J n{2e t )J [e t ) sin 2 6 t e t dt (44) 



are the transition probabilities, and 

oo oo 

" [2n(0) + l]J(fl)fidn 





$ = / e t eft sin 2 t 



(45) 
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is the additional phase picked up during the sweep. This 
phase shift is due to the renormalization of the gap A 
due to the interaction between the qubit and the envi- 
ronment. 

Comparison of Eqs. (|41|l and (|23[) shows that the phe- 
nomenological Eq. ill'MI) is correct in the perturbative 
limit. Also, since wi and W2 contain l + n(2et) and n(2e t ) 
respectively, in the low-temperature limit (T rcs <C A) the 
decay rate w\ is finite, while the excitation rate u>2 is 
thermally assisted: W2 ~ exp(— 2A/T res ). 
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